THE VLASOV-POISSON-BOLTZMANN SYSTEM FOR SOFT 

POTENTIALS 



RENJUN DUAN, TONG YANG, AND HUIJIANG ZHAO 

Abstract. An important physical model describing the dynamics of dilute 
weakly ionized plasmas in the collisional kinetic theory is the Vlasov-Poisson- 
Boltzmann system for which the plasma responds strongly to the self-consistent 
electrostatic force. This paper is concerned with the electron dynamics of 
kinetic plasmas in the whole space when the positive charged ion flow provides 
a spatially uniform background. We establish the global existence and optimal 
convergence rates of solutions near a global Maxwcllian to the Cauchy problem 
on the Vlasov-Poisson-Boltzmann system for angular cutoff soft potentials with 
—2 < 7 < 0. The main idea is to introduce a time dependent weight function 
in the velocity variable to capture the singularity of the cross-section at zero 
relative velocity. 
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1. Introduction 

1.1. Problem. The Vlasov-Poisson-Boltzmann (called VPB in the sequel for sim- 
plicity) system is a physical model describing the motion of dilute weakly ionized 
plasmas (e.g., electrons and ions in the case of two-species) in the collisional kinetic 
theory, where plasmas respond strongly to the self-consistent electrostatic force. In 
physics, the ion mass is usually much larger than the electron mass so that the 
electrons move much faster than the ions. Thus, the ions are often described by 
a fixed background n\,(x) and only the electrons move rapidly. In this simplified 
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case, the VPB system takes the form of 

(1-1) ft/ + £-Vx/ + Vx0-V«/ = Q(/ ) /) ) 

(1.2) A x <f> = fd£-l, 4>{x) ->■ as \x\ ->• oc, 

(1-3) f(0,x,0 = fo(x,0- 

Here the unknown / = f(t,x,£) > is the density distribution function of the 
particles located at x = (xi, X2,xs) £ R 3 with velocity £ = (£i, £2, £3) £ K 3 at time 
t > 0. The potential function </> = (f>(t, x) generating the self-consistent electric field 
V x 4> in (jl.ip is coupled with /(i, x,£) through the Poisson equation (|1.2[) . where 
rib(x) = 1 is chosen to be a unit constant. 

The bilinear collision operator Q acting only on the velocity variable, cf. [4j [14] 
135], is defined by 

(1.4) Q(f,g)= [f \^-^qo(9)(f(Qg(0-m*)9(0)<^^, 

where (£,£*) and (£',£*), denoting velocities of two particles before and after their 
collisions respectively, satisfy 

£' = £-[(£-£.) •£"]£«;, £ = £. + [(£-£.)• we§ 2 , 
by the conservation of momentum and energy 

= ici 2 + ie*i 2 = ie'i 2 + iei 2 . 

Note that the identity |^- ^| = - ^| holds. 

The function |C~ C*l 7 9o(^) in dEH is the cross-section depending only on — 
and cos9 = (£ — • — and it satisfies — 3 < 7 < 1 and it is assumed to 
satisfy the Grad's angular cutoff assumption < qo(9) < C\ cos0|. The exponent 
7 is determined by the potential of intcrmolccular forces, which is called the soft 
potential when —3 < 7 < 0, the Maxwell model when 7 = and the hard potential 
when < 7 < 1 including the hard sphere model 7 = 1, qo(Q) = cont.| cos0|. For 
the soft potential, the case — 2 < 7 < is called the moderately soft potential 
and —3 < 7 < —2 very soft potential, cf. [35] . Notice that the Coulomb potential 
coincides with the limit at 7 = — 3 for which the Boltzmann collision operator 
should be replaced by the Landau operator under the grazing collisions (23] [31] [35] . 

It is an interesting problem to consider the nonlinear stability and convergence 
rates of a spatially homogeneous steady state M for the Cauchy problem (jl.ll) . 
(|1.2|) . (|1.3[) when initial data fo is sufficiently close to M in a certain sense, where 

M = (27r)" 3 / 2 e-l«l 2 / 2 

is a normalized global Maxwcllian. This problem was first solved by Guo [22] for the 
VPB system on torus in the case of the hard sphere model. Since then, the case of 
non hard sphere models, i.e. —3 < 7 < 1, under the angular cutoff assumption, has 
remained open. In a recent work [13j for the VPB system over the whole space, we 
studied the problem for the case < 7 < 1 that includes both the Maxwell model 
and general hard potentials. However, the method used in [13] can not be applied 
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to the soft potentials, mainly because the collision frequency ~ (1 + [£|) 7 with 
7 < is degenerate in the large-velocity domain. 

In this paper we shall study the problem in the soft potential case by further 
developing the approach in |13j with the following extra ingredients to overcome 
the specific mathematical difficulties in dealing with soft potentials: 

• a new time- velocity-dependent weight in the form of exp{(£) 2 [g + A/(l + 
t) ]} to capture the dissipation for controlling the velocity growth in the 
nonlinear term for non hard sphere models, and 

• a time-frequency/time weighted method to overcome the large- velocity de- 
generacy in the energy dissipation. 

The approach and techniques that this paper together with |13j developed can be 
applied to some other collisional kinetic models for the non hard sphere interaction 
potential when an external forcing is present. 

1.2. Main results. To this end, as in [22], set the perturbation u = u(t,x,£) by 

Then, the Cauchy problem (|l.ll) - ()1.3[) of the VPB system is reformulated as 

(1.5) d t u + £ • \7 x u + V x <j> ■ V ? m - -£ • V x (j>u - V x <j> ■ £M 1/2 = Lit + T(u, u), 

(1.6) A x (j)= / M 1/2 ud£, 4>{x) — > as |x| — > oo, 

(1.7) u(0,s,0 = u (x,0 = M-^ifo - M), 

where the linearized collision operator L and the quadratic nonlinear term T are 
defined by 

Lu = M"5 |g (m, M 1/2 uj + Q (m 1/2 u, Mj } , 
T(u,u) = M~iQ (M 1/2 u,M 1/2 u) , 
respectively. Here, notice that due to (jl.6|) . <p can be determined in terms of u by 

(1.8) 0(t, x) = -j±-r * x [ M^uit, x, dt 

A-k\x\ J R 3 

By plugging the above formula into the equation (|1.5|) for the reformulated VPB 
system, one has a single evolution equation for the perturbation u, cf. [12( [9], 

Before stating our main result, we first introduce the following mixed time- 
velocity weight function 

(1.9) «V, g (*,6 = (0 7T e <?>2 [ 9+ ^> 

where tgR, < q < 1, < A < 1, 0< < 1/4, and (£_} = (1 + |£| 2 ) 1/2 . Notice 
that even though w r _ q (t,£) depends also on parameters A and we skip them for 
notational simplicity, and in many places we also use q(t) to denote q + (1 + t ^ for 
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brevity. For given u(t,x,£) with the corresponding function <p(t,x) given by (|1.8jl . 
we then define a temporal energy norm 

(1.10) HW*)= J2 ll%»i-*,a(*>0^«W|| + llv^(t)||ir*, 

\a\ + \l3\<N 

where N > is an integer, and I > N is a constant. 

The main result of this paper is stated as follows. Some more notations will be 
explained at the end of this section. 

Theorem 1.1. Let -2 < 7 < 0, N > 8, £ > §, £ > 1 + max (n, f - ±j and 
q(t) = q + (xq^p- > with < g <C 1, < A C 1, and < ■& < 1/4. Assume that 
f = M + M 1 / 2 ^ > and 

(1.11) // M 1 / 2 m dxd£ = 0. 

There are constants eo > 0, Co > such that if 

(1.12) £ \\w w - i>q (0,0d%u \\ + 

\a\+\P\<N 

then the Cauchy problem (|1.5p . (|1.6|) . (|1-7|) 0/ £/ie VP5 system admits a unique 
global solution u(t,x,t;) satisfying f(t,x,£) = M + M 1 / 2 u(f, x, ^) > and 

(1.13) sup{|u||iv Ag (t) + (1 + *) S |«|jv^-i, a (*) 

t>0 l 

+ (l + t)t ||V^(t)|| ffW _i} < C e . 

Remark 1.1. T/ie condition (|1.11[) implies that the ionized plasma system with 
electrons and ions are neutral at initial time. Due to the conservation of mass for 
(|1.1[) or (|1.5p . the neutral condition (jl.lip holds for all time t > 0. Notice that 
p. lip together with ||(1 + |a;|)uo||zi < 00 can oe used to remove the singularity of 
the Poisson kernel and to induce the same time-decay rate as in the case of the 
Boltzmann equation; see also Remark lS. 11 Moreover, in (|1.12p . we put the extra 
space-velocity weight 1 + \x\ + |£| 2~ on «o in Z\-norm. This is necessary in the 
proof of Lemma \6.1\ concerning the desired uniform-in-time a priori estimate. 

Remark 1.2. The problem of the global existence for the very soft potential case 
—3 < 7 < —2 is still left open. Despite this, we shall discuss in the next subsections 
the main difficulty. 

We conclude this subsection by pointing out some possible extensions of Theorem 
11.11 First, the arguments used in this paper can be adopted straightforwardly to 
deal with either the VPB system with soft potentials on torus with additional 
conservation laws as in or the two-species VPB system with soft potentials as 
in [21] [37]. Next, it could be quite interesting to apply the current approach as well 
as techniques in [20j [31] to consider the Vlasov-Poisson-Landau system [23] , where 
Q in (jl.ip is replaced by the Landau operator 

Q L (f, /) = v e • U B L (^ - £*) [/(£*) Vf/(£) - /(£)Vf/(£*)] C,} , 
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where B L (£) is a non- negative matrix given by 

^•(o = (««-|^)ier +2 , 7>-3. 

Note that the soft potential for the Landau operator corresponds to the case — 3 < 
7 < -2. 

1.3. Strategy of the proof. We first recall several elementary properties of the 
linearized operator L. First of all, L can be written as L = — v + K , where 

i/=i/(o~(i+i£ir 

denotes the collision frequency, and K is a velocity integral operator with a real 
symmetric integral kernel The explicit representation of v and K will be 

given in the next section. It is known that L is non-positive, the null space of L is 
given by 

M = span {M 1 / 2 , iM 1 ' 2 (1 < i < 3), I^M 1 / 2 } , 
and L is coercive in the sense that there is a constant ko > such that, cf. [4lH8l 

(1.14) -/ uLudt>K [ K£)|{I-PH 2 d£ 

JR 3 Jl 3 

holds for u = where I means the identity operator and P denotes the orthog- 
onal projection from L 2 to AT. As in [THIHD], for any given any u(t, x, £), one can 
write 

' Pu = {a(t, x) + b(t, x) ■ f + c(t, a:) (|£| 2 - 3) } M 1 / 2 , 
o= / M 1/2 urf^= / M 1/2 Pud£, 

jE 3 JK 3 

/ ^M^tidf = / ^M^Pttd^, 1 < i < 3, 

Js? Jl 3 

c=s/ (lf| a -3)M 1 /a ud e=i/ (l^-ajMVapude. 

V JR 3 JR 3 

so that we have the macro-micro decomposition introduced in |19j 

(1.15) u(t, x, = Pu(t, .t, + {I - P}u(t, x, 0- 

Here, Pu is called the macroscopic component of u(t, x, £) and {I — P}u the micro- 
scopic component of u(t, x, £). For later use, one can rewrite P as 

{Pu = PqU © PlU, 
P u = a(t,x)M 1/2 , 
Piu = {6(t, ar) ■ f + c(t, s) (|£| 2 - 3) } M 1 / 2 , 

where Po and Pi are the projectors corresponding to the hyperbolic and parabolic 
parts of the macroscopic component, respectively, cf. [9]. 

To prove Theorcm ll.il wc introduce the equivalent temporal energy functional 

(1-16) £W*) ~ MW*). 
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and the corresponding dissipation rate functional 

(i.i7) v NAq (t)= J2 |^ 1/2 ^i-^(*'0^{i-PH*) 

\a\ + \0\<N 

]T ||<eHii8i-< >g (*,ow-p}«(*)| 



(1 + 



i«ii 2 + E liv,a«( a ,6, c )" 2 

|a|<JV-l 



We also introduce the time-weighted temporal sup-energy 

3 

(1.18) X NAq (t) = sup £ NAq (s) + sup (1 + s)^£ N ,e-i, q (s) 

0<s<t 0<s<t 

+ sup (l + s)t ||V^(s)|| 2 ^ JV _ 1 . 

0<s<* 

Here, the exact definitions of £N,i,q(t) is given by (|5.16[) in the proof of Lemma HT21 
The strategy to prove Theorem 11.11 is to obtain the uniform-in-time estimates 
under the a priori assumption that X^,i,q{t) is small enough over < t < T for 
any given T > 0. Indeed, 

• one can deduce that there is a temporal energy functional £N,e, q (t) satisfy- 
ing (|1.16[) such that 

(1.19) Jt £NAq ^ + rVn ^i^> - °' 

for < t < T, where T>N,e, q (f) is given by (|1.17p . and moreover, 

• by combining the time-decay property of the linearized system, one can 
prove that 

(1-20) XnaM < C {e 2 NAq + X NAq (t) 2 } 

for < t < T, where £N,e,q depends only on initial data uq- 
Therefore, by the local existence of solutions as well as the continuity argument, 
Xjsr,e,q{t) is bounded uniformly in all time t > as long as £N,£,q is sufficiently 
small, which then implies Theorcm ll.il 

Wc remark that this kind of strategy used here was first developed in [TT] for the 
study of the time-decay property of the linearized Boltzmann equation with external 
forces and was later revisited in |10j for further generalization of the approach. 

1.4. Difficulties and ideas. We now outline detailed ideas to carry out the above 
strategy with emphasize on the specific mathematical difficulties. Before that, it 
is worth to pointing out that for the VPB system with soft potentials, the only 
result available so far is [24] in which global classical solutions near vacuum are 
constructed. One may expect that the work [32] and its recent improvement [53] 
by using the spectral analysis and the contraction mapping principle can be adapted 
to deal with this problem. However, we note that when the self-induced potential 
force is taken into account, even for the hard-sphere interaction, with the spec- 
tral property obtained in [16], the spectral theory corresponding to [32] has not 
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been known so far, partially because the Poisson equation produces an additional 
nonlocal term with singular kernels. 

Fortunately, the energy method recently developed in [201 [H] works well in the 
presence of the self-induced electric field [22] [38l H2] or even electromagnetic field 
[2T| [30J [7]. But due to the appearance of the nonlinear term £ ■ V x (f>{I — P}it, the 
direct application of the coercive estimate (|1.14j) of the linearized collision operator 
L works only for the hard-sphere interaction with 7 = 1. Indeed with a bit weaker 
(softer) than hard-sphere interaction, such a term is beyond the control by cither the 
usual energy or the energy dissipation rate so that even the local-in-time solutions 
can not be constructed within this framework. 

To overcome this difficulty, a new weighed energy method was introduced in |13j 
by the authors of this paper to deal with the VPB system with hard potentials. 
Such a method is based on the use of a mixed time-velocity weight function 

(1.21) ^(U) = (0 4 e^, 

where £ E M, A > and > are suitably chosen constants. One of the most 
important ingredients is to combine the time-decay of solutions with the usual 
weighted energy inequalities in order to obtain the uniform-in-time estimates. 

The main purpose of this paper is to generalize the above argument so that it can 
be adapted to deal with the soft potential case. Compared with the hard potential 
case, the main difference lies in the fact that at high velocity the dissipation is much 
weaker than the instant energy. Our main ideas to overcome this are explained as 
follows. 

Firstly, similar to the hard potential case [13] , we introduce an exponential weight 
factor 

< 9 (i,£) = exp j<0 2 

in the weight function w T , q {t,£,) given by (jl.9j) . Notice that q(t) := q + A/(l + t) 
satisfies < q(t) <C 1 by choices of q, A and ?9 as stated in Theorem 1 1.1[ and hence 
all the estimates on v and K obtained in (31] can still be applied with respect to the 
weight function w T>q (t, £) considered here. A key observation to use the exponential 
factor is based on the identity 

w e T>q (t,S)d t g(t,x,0 = d t {w e T>q (t,Og(t,x,0} + ™® w e Ttq (t,Og(t,x,®. 

Thus one can gain from the second term on the right hand side of the above identity 
an additional second-order velocity moment (£) 2 with a compensation that the 
magnitude of the additional dissipation term decays in time with a rate (l+t)~( 1+ ^ . 
With this, the nonlinear term £■ \7 x (fi{I— P}u can be controlled as long as the electric 
field V x 4> has the time-decay rate not slower than (1 +t)~( 1+i> \ Notice that in the 
case of the whole space, || V x 4>\\ h n~i as a part of the high-order energy functional 
decays at most as (1 + and thus it is natural to require < i? < 1/4. 

Moreover, as used for hard potentials in |13j . the first-order velocity moment (£) 
in the exponential factor of (|1.21j) is actually enough to deal with the large- velocity 
growth in £ • V x <f){l — P}u. The reason why we have chosen (£) 2 for soft potentials 
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is that one has to control another nonlinear term V x 4>- V^{I — P}u. We shall clarify 
this point in more details later. 

The second idea is to control the velocity £ derivative of u. It is well-known 
that Vjm may grow in time. To overcome such a difficulty, we apply the algebraic 
velocity weight factor, introduced in [20 ] [T8 ] . 

= (0^ m - e \ e>\P\, 

in the weight w Ti q(t,£) given in (|I.9[) with r = \/3\ —£. This algebraic factor implies 
two properties of the velocity weight: one is that wfjgi^O > 1 holds true due to 
t > an d the other one is that the higher the order of the ^-derivatives, the more 
negative velocity weights. The restriction £ > |/3|, more precisely 

^-l>max|iV,|-i|>|/3|, 

results from the fact that one has to use the positive-power algebraic velocity weight 
w ?p\—t(Q so as to obtain the closed estimate (|1.20[) on Xjv/, g (t), because not only 
the nonlinear term £ ■ V x <fiu contains the first-order velocity growth but also initial 
data is supposed to have an extra positive-power velocity weig ht in the 

time-decay estimate (|3.5j) for the evolution operator of the linearized VPB system. 

As pointed out in [18] for the study of the Boltzmann equation with soft poten- 
tials, the dependence of the weight wj^,^) on |/3| makes it possible to control the 
linear term £ • V x u in terms of the purely x-derivative dissipation for the weighted 
energy estimate on the mixed x-£ derivatives. However, the algebraic velocity factor 
w \p\-l(£) produces an additional difficulty on the nonlinear term V ' x <j> ■ V^it in the 
presence of the self-consistent electric field for the VPB system with soft potentials. 
To obtain the velocity weighted derivative estimate on such a nonlinear term, one 
should put an extra negative-power function (£) 7 in front of V^{I — P}u in the 
nonlinear term V ' x <f> ■ V^u so that the velocity growth (£) -7 comes up to have a 
balance. Then, as long as 

-2 < 7 < 0, 

it is fortunate that the dissipation functional D^^^it) given by (|1.17[) containing 
the second-order moment (£) 2 can be used to control the term V ' x <j> ' V^u. 

At this point, it is not known how to use the argument of this paper to deal with 
the very soft potential case —3 < 7 < —2. However, inspired by [23], we remark 
that one possible way is to have a stronger dissipation property of the linearized 
collision operator L, particularly the possible velocity diffusion effect. In fact, 
for the Boltzmann equation without angular cutoff assumption, from the coercive 
estimate on the linearized Boltzmann collision operator L obtained in [T] [T7], it 
seems hopeful to deal with the case —3 < 7 < —2 for such a physical situation, 
which will be reported in the future work. 

Another ingredient of our analysis is the decay of solutions to the VPB system 
for soft potentials. Recall that the large-time behavior of global solutions has been 
studied extensively in recent years by using different approaches. One approach 
which usually leads to slower time-decay than in the linearized level is used in [38] 
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on the basis of the improved energy estimates together with functional inequali- 
ties. The method of thirteen moments and compensation functions proposed by 
Kawashima in [27] which gives the optimal time rate without using the spectral 
theory; see |15j and |37j for two applications. Recently, concerning with the opti- 
mal time rate, a time- frequency analysis method has been developed in (9J [8] Ej- 
Precisely, in the same spirit of [36], some time- frequency functionals or interactive 
energy functionals are constructed in [9] E] to capture the dissipation of the 
degenerate components of the full system. 

Back to the nonlinear VPB system, the main difficulty of deducing the decay 
rates of solutions for the soft potential is caused by the lack of a spectral gap for the 
linearized collision operator L. Unlike the periodic domain [31 j . we need a careful 
and delicate estimate on the time decay of solutions to the corresponding linearized 
equation in the case of the whole space K 3 . Following the recent work [5] [5], 
our analysis is based on the weighted energy estimates, a time-frequency analysis 
method, and the construction of some interactive energy functionals, which gives a 
new and concise proof of the decay of the solution to the linearized equation with 
soft potentials. Here, we should mention the recent work [29j by Strain. Different 
from this work, by starting from the inequality 

dtEt(u) + Kp(k)Et-i(u) < 0, 

we use a new time-frequency weighted approach together with the time-frequency 
splitting technique in order to deduce the time-decay estimate on the solution u for 
the linearized VPB system; see Step 3 in the proof of Theorem 13 . 1 1 and the identity 

Finally we also point out a difference between the soft potential case and hard 
potentials studied in |13| for obtaining the time-decay of the energy functional 
^N,e,q{t)- The starting point is the Lyapunov inequality (| 1 . . In the hard poten- 
tial case, the Gronwall inequality can be directly used when the lowest-order terms 
|| (a, 6, c)|| 2 and ||V K 0|| 2 are added into T>^^^ q (t). However, this fails for the soft 
potentials because T>N,t, q {t) given by (|1.17[) contains the extra factor <~ (£) 7 
which is degenerate at large-velocity domain when 7 < 0. To overcome this dif- 
ficulty, we have used the time-weighted estimate on (| 1 . 19|) as well as an iterative 
technique on the basis of the inequality 

^v,i, 9 (*) + ll(MWII 2 + l|v*#)|| 2 > *£ N ,i- U (t) 

for any £; see Step 2 in the proof of Lemma [6. II for details. 

Before concluding this introduction, we mention some previous results concern- 
ing the study of the VPB system in other respects, cf. [3] [35], the Boltzmann 
equation with soft potentials, cf. [33J[5][5], and also the exponential rate for the 
Boltzmann equation with general potentials in the collision kernel but under addi- 
tional conditions on the regularity of solutions, cf. [6]. 

The rest of this paper is organized as follows. In Section [21 we will list some 
estimates on v and K proved in [31] . In Section [3] we will study the time-decay 
property of the linearized VPB system without any source. In Section [4] we give 
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the estimates on the nonlinear terms with respect to the weight function w T , q (t, £). 
In Section[5l we will close the a priori estimates on the solution to derive the desired 
inequality (|1.19[) . In the last section, we will prove (|1.20[) to conclude Theorcm ll.il 

1.5. Notations. Throughout this paper, C denotes some positive (generally large) 
constant and k denotes some positive (generally small) constant, where both C and 
k may take different values in different places. A < B means A < and A ~ B 
means kA < B < —A, both for a generic constant < K < 1. For an integer m > 0, 
we use H™ £ , H™, to denote the usual Hilbert spaces H m (B? x x H m (R*), 
H m (M.^), respectively, and L 2 , L 2 , L 2 are used for the case when m = 0. When 
without any confusion, we use H m to denote H™ and use L 2 to denote L 2 or L 2 X ^. 
We denote (■, ■) by the inner product over L 2 For q > 1, we also define the mixed 
velocity-space Lebesgue space Z q = L 2 (L q x ) = L 2 (R|; L q (R x )) with the norm 

Hk=(7 (I Wx,0\ q dx) /q dti) , u = u(x,£)eZ q . 



For multi-indices a = (ai, 02,03) and (3 = (/3 1; /3 2 , A3), we denote dp = d x d^, that 
is, <9^ = d^d^d^d^d^d^. The length of a is \a\ = a x +a 2 + a 3 and the length 
of k \P\=Pi+p 2 +pl 



2. Preliminaries 
Recall that L = — v + K is defined by 

K0= // |e-^| 7 9oWM(^)da;de*~(l + iei) 7 , 



and 



x-u(0 = // ie-e.r 8 d(fi)M i /2(e,)M i /2(^)„(f)d W de, 

J Jl 3 xS 2 
J JR 3 x S 2 

- ff \^-^q (e)M 1 / 2 (^)M 1 / 2 (Ou(^)du J ^ 

J JR3 x S 2 



We list in the following lemma velocity weighted estimates on the collision frequency 
j^(£) and the integral operator K with respect to the velocity weight function 

u>r,q(0 = (0 7T e q{ ° 2 , t£R, < q < 1. 
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Lemma 2.1 ([31]). Let -3 < 7 < 0, r G K } and < g < 1. 7/ |/3| > 0, f/ien /or 
any 77 > 0, there is C n > such that 



17 



\P1\<\P 

If \P\ > 7 i/ien /or any r\ > 0, i/iere is > smc/j i/iai 



(2.1) 



t4 tS (£)dp(Kf)gdt 



< 



-a 



xi ? i<2c,(0' 7T i/r^ 



E 

l/3i|<l/9| 
2 7 r| f |2 ( 



u{0wl^)\d Pl f\ 2 di 



KC)<,-(0lffl 2 ^ 



For later use, let us write down the macroscopic equations of the VPB system 
up to third-order moments by applying the macro-micro decomposition (|1.15[) in- 
troduced in [19]. For that, as in [9], define moment functions Oy(-) and A^(-), 
1 < i,j < 3, by 

(2.2) Qij(v)= [ {^-l)M^ 2 vd^, Ai(v) = ±-[ m 2 -5)^M 1/2 vd^ 

for any u = v(£)- Then, one can derive from (|1.5|) - (|1.6j) a fluid- type system of 
equations 

( d t a + V x ■ b = 0, 



(2.3) 



and 



(2.4) { 



d t b + y x (a + 2c) + V* • 6({I - P}n) - V x = V X <M, 

ftc + ^V* • b + V x • A({I - P}u) = ^V x <j> ■ b, 
, = a, 
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9 t e y ({i - p} u ) + ^ + - -s, u v x ■ b - — • a({i - p} u ) 



%(r + G) - -SijV^-b, 



{ d t Ai({I - P}u) + d lC = A,(r + G) 



with 



r = -£ ■ V X {I - P}n + Lu, G = T(w, «) + -£■ V x </>u - V x • V^it, 

where r is a linear term related only to the micro component {I — P}n and G 
is a quadratic nonlinear term. Here and hereafter, for simplicity, we used dj to 
denote d Xj for each j = 1, 2, 3. The above fluid-type system (|2.3|) - (|2.4[) plays a key 
role in the analysis of the zero-order energy estimate and the dissipation of the 
macroscopic component (a, b, c) in the case of the whole space; see [El [El El U\- 
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3. Time decay for the evolution operator 
Consider the linearized homogeneous equation 

(3.1) d t u + £ • V x u - V x <j) ■ £M 1/2 = Lu, 
with 

(3.2) = ^—-*x( M 1/2 u(t,x,£_)d^ ->• as \x\ -> oo. 

4?r|a;| J R3 

Given u — uo(x,£), e tB u is the solution to the Cauchy problem (|3.1[) - (|3.2jl with 
u\t=o = uq. For an integer m, set the index <j m of the time-decay rate by 

3 m 
° m = 4 + 7' 

which corresponds to the one for the case of the usual heat kernel in three dimen- 
sions. 

The main result of this section is stated as follows. 

Theorem 3.1. Set \i = £t(£) := Let -3 < 7 < 0, ^ > 0, a > 0, m = \a\, 

£q > 2(7 m , and assume 

(3.3) / cic,dx = 0, / (1 + |a;|)|ao| dx < 00, 

JR3 J R 3 



tin 



d 



(3.4) \\n e+e °uo\\ Zi + \\^ +to d a u \\ < 00. 
Then, the evolution operator e tB satisfies 

(3.5) \\^d a e tB u \\ + H^VaA^Poe^uoll 

< C7(l + 1)-°™ {]\^ +t °u4 Zi + W^'d-uoW + ||(1 + \x\)ao\\ Li ) 

for any t > 0. 

Proof. The proof is divided by three steps. 

Step 1. The Fourier transform of (|3.1l) - (|3.2p gives 

(3.6) d t u + i£ ■ ku - ik(f> • £M 1/2 = Lu, 

with — \k\ 2 (j> = a. By taking the complex inner product of the above equation with 



u, integrating it over R| and using (JTTT4J) , as in [9], one has 



< 0. 



(3-7) d t nu\\i l + L^j+K ^{I-P} 

Furthermore, similar to derive (|2.3p and (|2.4p from the nonlinear VPB system, one 
can also derive from (|3.1[) - (|3.2[) a fluid- type system of linear equations 

d t a + V x ■ b = 0, 

d t b + V x (a + 2c) + V, • 9({I - P}u) - V x </> = 0, 
d t c + i V, • b + ^V x ■ A({I - P}u) = 0, 

A x <f> = a, 
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and 



d t & t] ({I - P}u) + dibj + djbi - -%V X • b 



10 



-S ij V x -A({I-P}u) = Q ij (r), 

d t A i ({I-P}u) + d i c = A i (r) 

with r = — £ • V X {1 — P}u + Lu, where O(-) and A(-) are defined in (|2.2[) . From the 
above fluid- type system, as in [TJ, one can deduce 



(3.8) d t $l£ int (ii{t,k)) + k 
where £ lnt (u(t, k)) is given by 



[j^L \M~c)\ 2 + |a| 2 } < C \v^{l P\z 



£ int (u(t,k)) 



1 



(ike | A({I-P}u)) 



^ [ikjb m + ik m bj - ^5 jm ik ■ b \ Q jm {{! - P}u) 

jm 



+Ki ( ika I b 



Therefore, for < K2 <C 1, a suitable linear combination of (|3.7j) and (J378J) gives 

IAI2 



(3.9) d t 



\k\- 



+ K 2 ^£ m \u{t,k)) 



+ K 



*// 2 {I-P}? 



a i + |fc| 5 



[a,b,c) 



a\ z y < 0. 



Here, notice that since 

\£ int (u(t,k))\<C{\\{I-P}u\\U 



(a, &, c) 



we have taken n 2 > small enough such that 

|2 



(3.10) 



u\\U + tt^+ n 2 ^£ m \u(t, k)) ~ \\u\\i 2 + 



\k\> 



lL 'i 1 |fc| 2 

Step 2. Applying {I - P} to (|3^]l . we have 

<9 t {I - P}u + if; • fc{I - P}u = L{I - P}ii + Pi£ • M - if • kPu. 

By further taking the complex inner product of the above equation with /z 2£ (£){I 
P}u and integrating it over R|, one has 

(3.ii) \d t - p}«|| a + k - P}«||i2 

< / (K{\ - P}u | M 2 ^(e){I - P}u) d£ 

+ 3i [ (Pif • M - if ■ kPu | /r"(f){I - P}w) eZf. 

JR 3 
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It follows from (|2.1[) that for any 77 > 0, 



(3.12) 



/ (K{I - P}u I M 2 ^){I - P}u) d$ 

Jr 3 

< jr?p/y{i-p} 



X|c|<ac,{0"^{I-P}fi|} 

vV{i-p}fi 



which further by the Cauchy-Schwarz inequality implies 



/ (K{I P}u I ^(C){I - P}u) d£ 

JR 3 



< 



i^/V{i-PK 



^{I-Pjf 



Notice that over fc| < 1, 



/ (Pi£ • few - if • fcP« I - p }«) d?X|*|<i 

^/ 2 {I-P} 



< c 



(a, b, c) 



1 + 1*1 



Plugging the above two inequalities into p. lip and fixing a small enough constant 
77 > give 



(3.13) d t ||//{I - P}w|| i2 X |fc|<i + « " P}«|| L 2 X|fc|<i 



< C 



/^{I-P}* 



l\ l + \k\ 



k\ 2 — , 2 
(a, b, c 



To obtain the velocity-weighted estimate for the pointwise time-frequency vari- 
ables over |fc| > 1, we directly take the complex inner product of (|3.6p with /j, 2e (^)u 
and integrate it over R| to get that 



1 



(3.14) 1^11^11' + i/Vy^.p^ 

= -» / (i/(C){I-P}u I A» M (0P«) + » / (K{I-P}u\fi 2e (Qu) <% 

Jm Jr 3 

+ U I (ik<j> ■ (M 1 ' 2 I ^ 2£ {Ou) d£. 

Here, Cauchy-Schwarz's inequality implies 

-u f - p}^ 1 m 2 '(0p«) # < c ( v^ii - p} ? 

Jr 3 \ 



Li 



(a, b, c) 
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and 



M 3 



= SR / (ifc0 • CM 1 / 2 | //"(£) [Pu + {I - P}u]) d£ 



< c 



v 1/2 {l-V}u 



(a, &, c) 



where |fc| 2 |</>| 2 = |a| 2 /|/c| 2 due to the Poisson equation — \k\ 2 cj) — a was used. From 
(12. ip , similar to the proof of (|3.12[) , one has 

/ (K{I-P}u\^(Ou) dd 

[ (A-{I-P}fi| M M (Opfi + {I-P}fi]) 

JR 3 



< 



/{I-P}i 



,1/2 



{I-P}t 



(a, 6, c) 



By plugging the above two estimates into (|3.14j) . fixing a small constant r] > and 



2 1 k i ~ 

using 1 < pjjrcp X|fc|>i; it follows that 

(3.15) 9 t ||m £ ^||l2 X|fc|>i + K H^ -1 ! 1 ~ p I"||l 2 X|fe|>i 



< C 



^{I - P} t 



1 + |*| 



fc| 2 , ! 

(a, o, c 



Now, for properly chosen constants < K3, K4 -C 1, a suitable linear combination 
of (HU), (f37T3|) and (f3TT5]) yields that whenever i > 0, 

(3.16) d t E e (u) + nD e {u) < 0, 

holds true for any t > 0, fc G R 3 , where Eg{u) and Di{u) are given by 

£i(tt) = |H|i 2 + |^L + /c 2 Kr nt (u(t,fe)) 

5 |fc| 

+ K 3 - P}u|L 2 X|fe|<l + K4 I|/«||r2 X|fc|>l, 

|fc| 2 



(o,6,c) + |a|" 



lL I 1 + |fc| 2 

Due to (|3.10p and the fact that Pit decays exponentially in £, it is straightforward 
to check that for t > 0, 



jar 



Step 3. To the end, set 



p(k) 



\k\< 



l + \k\ 2 ' 

Let < e < 1 and J > be chosen later. Multiplying ([3~16f by [1 + ep(k)t] J , 

(3.17) ^{[l + ep(k)t} J E e (u)} + n[l + ep{k)t} J D e {u) < J[l + ep(fc)t]' / - 1 e i o(fc)^(u). 
at 
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To estimate the right-hand term, we bound Ee(u) in the way that 



E t (u) ~ \\n u\\ L , + 

n f-n 2 / \ |a| 2 

= ||M U\\ L 2 (X|fc|<l+X|fe|>l) + TjT2 

< (||^{! - P HIl 2 + ||//Pu|| i2 )X|fc|<l + H^ILa X|fe|>l 



jar 

2 



That is, 



< || A* {I - p }"|| L 2 X|fe|<i + ||m "|| L 2 X|fe|>i + | (a, fo, c) | + |— 2- 



^E[(u)+E n (u), 

E l (U) = ||/{I - P }"||i2 X|fc|<l + ||^W|| L 2 X|fe|>l, 



E n (u) = \(a,b,c)\ 2 + ^. 

Therefore, for the right-hand term of (|3.17p . one has 

J[l + ep(fc)i] J " 1 e / 9(fc)^(w) 
< J[l + ep{k)t] J - 1 ep{k)El{u) + J[l + ep(k)t} J - 1 ep(k)E II (u) := ife + R 2 

We estimate R\ and R 2 as follows. First, for R 2 , 

R 2 = J[l + ep{k)t] J - l ep{k)E u {u) 

= J[l + ep(fc)t] J - 1 ep(fc)( |(^)| 2 + }^L) 



= e J[l + ep(k)t] 



J-i 



|fc| 2 2 



l + |*| s 



(a, 6, c) 



< eJ[l + ep(k)t] J D e (u). 
Thus, one can let 



l + |fc| 2 

|2 "1 



£J< 4 

where k > is given in ([336]) . such that 

i? 2 < ^[l + ep(fc)t] J D*(«)- 

For we use the splitting 

( 3 - 18 ) 1 = X^(i)<[l+ep(k)t] + X M 2(£)>[i +ep (fc) t ], 

so that 

Ej(u) = E}(ux^(£)<[i +ep (k)t]) + E e(u Xl ,2 {()>[1+ep{k)t] ) := S/ < (tt) + 



with 
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E^iu) = I M 2 '(OI{I - P}«| 2 dtx\k\<i 

J^-(i)<[l+ep(k)t] 

+ / ^(0\u\ 2 d^ X \k\>i, 

J^(t)<[l+ep(k)t] 

E I l > (u)= f ^(0|{I-P}fi| 2 dfr|*|<i 

H 2l {Q\u\ 2 d£, X \k\>i- 



V 2 (?)>[i+ep(fc)t] 

Thus, with this splitting, R\ is bounded by 

Ri = J[l + ep{k)t] J - l ep{k)El{u) 

< J[l + e i o(fc)t] J - 1 ep(fc){^ /< (t/) + 

< eJ[l + ep(fc)t] J ' p{k)E I t < 1 {u) 

+J[1 + ep(k)t} J - 1 ep(k)E I e > (u) := R n + R 12 , 

where we have used 

El<{u) < [l + epQe^E^u). 

For Rn, notice that 

pifyE^iu) < p(fc)(||/- 1 {I - p }u\\ 2 L 2 X\k\<i + II^^H^ X|fe|>i) 

< H^-Hl-Plwll' + p{k) KM* 

< C x D t {u) 

for a generic constant C\ > 1, and hence 

i?n < C ie J[l + ep(k)t] J D e (u). 
One can make eJ further small such that 



and hence 



Cl eJ<^ 



Rn < - [l + ep(k)t] J D e (u). 



To estimate R12, let p > 1 be chosen later and compute 

R 12 = J[l + ep(k)t]-Pep(k) • [1 + epik^+P^E^iu) 

<J[l + ep{k)t]-Vep{k)-El>j+ 

Noticing that the estimate 

E e+j +p -i(ti) < ^/+j+p_i(w) < < ^+j+p-i("o) 

holds true by (|3.16p due to £ + J + p — 1 > 0, one has 

R12 < J[l + ep(/c)t]- p ep(/c)^ + , /+p _ 1 (u ). 
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Therefore, collecting all estimates above, it follows from (|3.17|) that 

D t (u) 

< J[l + ep(^^ p e^)£^fJ+P-lM• 



|{[l + ep(k)t] J E e (u)} + |[1 + ep(k)t] J 'D t (u) 



Integrating this inequality, using 

/ [1 + ep(k)s]- p ep(k)ds < (1 + = C p < oo 

Jo Jo 

for p > 1, and noticing J + p — 1 > 0, one has 

[1 + ep{k)t] J E e (u) < E e (u ) + JC p Ee+ J+p -i(u ) < (1 + JC p )E e+ j +p ^(u ), 
that is, whenever I > 0, 

£ [1 + ep(fc)t]- J ^+, /+p _i(wo), 
for any t > 0, fcgl 3 , where the parameters p, e and J with 
p > 1, < e < 1, J > 0, CieJ < |. 

are still to be chosen. 

Now, for any given £q > 2a m , we fix J > 2a m , p > 1 such that £q = J + p — 1 to 
have 

Ei{u) < [l + ep{k)t}~ J E Wo (u ). 
Since J > 2o~ m , in the completely same way as in [5] and |13j by considering the 
frequency integration over R| = {|fc| < 1} U {|fc > 1} with a little modification of 
the proof in [251 [27] , one can derive the desired time-decay property (|3.5p under 
conditions (|3.3|) and (|3.4[) ; the details are omitted for brevity. This completes the 
proof of Theorem 13.11 □ 



Remark 3.1. The rate index a m in (|3.5j) is optimal in the sense that it is the same 
as in the case of the Boltzmann equation. In fact, from |32[ I34j , one has 



^a e t(-?-v x +L; 



"o 



<c(i + t)-*i«i (Hlk + ll^iill), 



for any t > 0, where L the linearized Boltzmann operator in the case < 7 < 1. 
Here, different from the hard potential case, the extra velocity weight (£) ~ in 
(|3.5[) on initial data uq for soft potentials — 3 < 7 < is needed. Moreover, as 
we mentioned before, the condition (|3 - 3[) is used to remove the singularity of the 
Poisson kernel and recover the optimal rate index o~ m , otherwise, only the rate 
(1 + t)~i^ ~ can be obtained, see [5] and references therein. 

4. Estimate on nonlinear terms 

This section concerns some estimates on each nonlinear term in G = T(u, u) + 
' V ' x4>u ~ W X <P • V^w which will be needed in the next section. For this purpose, to 

simplify the notations, in the proof of all subsequent lemmas, w T is used to denote 

w r,q{t, £)• That is 

W T = W T , q (t,0 = <0 7r e 9(t) ^ 2 , q(t) = q + — 
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where < q(t) <SC 1 is fixed and r = \j3\ — i < depending on the order of velocity 
derivatives may take different values in different places. Moreover, we also use in 
the proof 

U\ = Pm, U2 = {I — P}u. 

Now we turn to state the estimates on the nonlinear terms. To make the presen- 
tation clear, we divide it into several subsections. The first subsection is concerned 
with the estimates on T(u,u). 

4.1. Estimate on T(u,u). We always use the decomposition 

(4.1) r(u, u) = r(Pu, Pu) + r(Pu, {i - p}u) + r({i - p}u, p u ) 

+ r({i-p}u,{i-p}u). 

Recall 

r(f,g)=T+(f,g)-T-(f,g) 

= [ le-^PM 1 / 2 ^) [/ qo(6)f(£)g(Odu 

JR 3 US 2 



5(0 / le-^rM 1 /^) [/ q^doj 

JR3 JS 2 



fit*) da* 



Our first result in this subsection is concerned with the estimate on it) both 
without and with the weight. 

Lemma 4.1. Let N > 4, £ > 0, < q(t) < 1. Jf ZioZds £/ia£ 



(4.2) (T(u,u),u) < C||(a,6,c)|| Hl \\V x (a,b,c) 



/Va{i_p} t 



+ c\\\y x (a,b 1 c)\\ H1 + ]T Pp^-pH} ^ 1/2 {i-Ph 

|a] + |/3|<4 



(4.3) (r(«, M ), w 2 _^(u){i-P} u ) 

<C||(a ) &,c)|| ff i ||V>,MII •||^/ 2 {I-P} 

+ C {E N , Lq {t)} 1/2 V V 2 w-t,Si 0{I - P}« 



Proof. We will prove only (|4.3j) . since (|4.2j) can be obtained in the similar way by 
noticing (T(u, u), u) = (T(u, u), {I — P}u). To this end, set 

h = (T(u,u),w 2 _ itq (t,0{l-'P}u) , 

and denote /i,n, Ii,i2, Ii,2i, -^1,22 to be the terms corresponding to the decomposi- 
tion (|4.ip . Now we turn to estimate these terms term by term. 

Estimate on Ii t u: Since T(ui,Ui) decays exponentially in £, 



/i,ii<C / |(a,6,c)| 2 

JR3 



dx < C\\(a,b,c)\\ H i \\V x (a,b,c) 
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where Holder's inequality with 1 = 1/3 + 1/6 + 1/2 and Sobolev's inequalities 

\\(a,b,c)\\ L » < C||(o,6,c)||hi, ||(o, b, c)\\ L , < C\\V x (a,b,c)\\ 
have been used. 

Estimate on I\ t \2'- For the loss term, 



^12 < C 



\u 2 (0\-\(a,b,c)nO-w 2 _ e (Z)\u 2 (Z)\dxdt 



< C sup \(a,b,c)\ ■ u^ 2 w-i(0u2 <C\\V x (a,b,c)\\ H i ^ 1/2 w_,(£)m2 

x£M 3 

where we have used in the first inequality 



(4.4) 



f ie-e.rM«'/ 2 K0de.<c<0' 

JM. 3 



for q' > and —3 < 7 < 0. For the gain term, using Holder's inequality, lti 2 1S 
bounded by 



(4.5) 



1^ - ^\ y q (e)M 1 ^(^)w 2 _ e (0\u 1 (OMO\ 2 dxdtdt.du 



x 1/2 



K - e*r9oWM i/2 (c*)^,(e)iw 2 (e)i 2 dx^ m duj 



s 1/2 



From (|4.4[) , the second term in (|4.5[) is bounded by ||^ 1 / 2 w_f(^)M2||. By noticing 
(0<min{(e:)(O,(e> + (O}and 

(4.6) «;_/(£) = (e)-T' e «W<« 

for £ > 0, the first term in (|4.5j) is bounded by 



< c 
= c 
= c 



\t-^qo(e)MH^)w^ e (m q \Q\(aAc)\ 2 \u 2 (0 



1/2 



ie-€*r<zo(^(e')M^(e)iK&, C )i 2 K(ni 2 



1/2 



1/2 



1/2 



\C-^Qo(0)wl e (OM^^)\(a,b,c)\ 2 \u 2 (O\ 2 

where we have used |Pu(£)l ^ C|(a, 6, c) | IVE 9 '/ 2 (^^ ) with < q' < 1 in the first 
equation, M 1 / 2 ^*) < C, the inequality gU) and W-i(^)M q ' ' /2 (C) < C in the 
second equation, the identity |£ — £»| = |£' — in the third equation and the 
change of variable (£, £) — > (£', with the unit Jacobian in the last equation. Thus, 
from (|4.4p . the first term in (|4.5[) is bounded by Csup^. \ {a, b, c)| ■ ||^ 1 / 2 -u;_f (£) u 2 • 
Therefore, 

It 12 < C\\V x (a,b,c)\\ H i v 1/2 w-t{0u2 
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Estimate on Ix t 2i- This can be done in the completely same way as for Ji,i2, so 
that 

h,2i <C\\V x (a,b,c)\\ m jiA' 2 w-i(Z)u2 

Estimate on l\.2i'- By [311 Lemma 3, page 305], 

(4.7) Ji,22 < C {£ NJ . q (t)} 1/2 |p /2 W _,{I - P}« 

and for brevity we skip the proof of this term. 

Now, the desired inequality (|4.3|) follows by collecting all the above estimates. 
This completes the proof of Lemma |4T1 □ 



For the weighted estimates on d^T(u, u), we have 

Lemma 4.2. Let N > 8, 1 < \a\ + |/3| < TV, £ > \/3\ and < q(t) < 1. For given 
u = u(t, x, £), define u a p as u a p = d a u if = and u a fj = dp{I — P}u if |/3| > 1. 
Then, it holds that 



v 1/2 u af3 



(4.8) ^T(u, U ),w 2 0l _ itq (t,Ou a p 

< C\\(a,b,c)\\ H N \\V x (a,b,c)\\ HN _ 1 
Proof. Take a, with 1 < \a\ + \(3\ < N. Write 



+ C{£ NAq (t)} 1/2 V NAq (t). 



where the summation is taken over /3q + fi\ + /?2 = ft and a\ + ct2 = a, and To is 
given by 



To (d£f,d%g) = T+ (dfj,d%g) - T (dfj,d%g) 



3 xS 2 

In what follows we set 

h = 



3 xS 2 



(d^u,d^uj ^^(i, 



and denote the terms corresponding to the decomposition (|4.1[) as 22. n, ^2.12, ^2,21 
and /2,22- These terms will be estimated term by term as follows. 

Estimate on Ja.ii* Notice that To (dp*u±, d^uij decays exponentially in £ since 
-3 < 7 < 0. Then, 



h,u<C / \d a Ha,b,c)\ ■ \d a *(a,b,c)\ 



v 1/2 u af} 



dx 



< C\\(a,b,c)\\ H N ||V x (a > 6,c)|| ffJ ,_i 



v 1/2 u aj3 
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Estimate on For the loss term, since 



-*2,12 ^ 



|£ - 6|>(0)M%) {M%)|d Ql (a, &, c)|} 
x d%u 2 (Z)\-wfa_ e (Z)\u a0 (O\ 
<C [dx \d a ^(a,b,c)\ ( d^(0 J w^_ e (0 0£u a (O«a/s(O 



we have for |«i| < N/2 that 



I 2 ~ 12 < C sup \d a ^a,b,c)\ • ^^w^iOd^uz 



v 1/2 w m _ e {C)u a p < C{£ NAq (t)} 1/2 V NAq (t), 



while for the case \ai\ > N/2, we can get that 



I- 12 <cjdx \d^{a,b,c)\ ■ jv^ww-tWdfiuaW^ \\v 1/2 w w _ e (Ou a p 
<Csup uWww^Qdgui 2 -\\d a H^b,c)\\- ^ 2 w m _ e (Ou af! 

< C {£ N , e , q (t)} 1/2 V N . Lq {t). 



For the gain term I 2 i2i 



^12 < / |£-^|\o(0)M%){M%)|^Ac)|} 

x ^>(£0|-^ M (OK/?(OI- 



We now turn to deduce an estimate on the right-hand side of the above inequality 
by considering the two cases |ai| < N/2 and |«i| > N/2, respectively. In the case 
of |ai| < N/2, we have from Holder's inequality that 



(4.9) I+ 12 < J dx \d ai (a,b,c) 



1/2 



1/2 
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As for estimating the first factor in (|4.5[) . we have 



< c 



which is bounded by 



Moreover, 



< C 



< c 



^ 1/2 U-| /3 |_ € (^)w ct ^ 



Therefore, plugging the above two estimates into (|4.9I) . in the case when |ai | < N/2, 
we have 

7+ 12 < Csu V \d^(a,b,c)\ ■ ^ 1/2 w w _ i (Od^u 2 

X 

v 1/2 w w _ e (Ou a p\\ < C{£ N . e . q (t)} 1/2 V N . Lq (t). 

The discussion for the case |ai| > N/2 is divided into the following three subcases: 
In Di = > ±|£|}, we have 



r+ <c I dx \d a Ha,b,c)\ 



Here, by Holder's inequality, we can deduce that 

\^-^qo(0)M^(OM^(^)\d^u 2 (C)u a p(O 

ie-e*r9oWM^(OM^(^)| u ^(oi 2 



< 



< c 



u 1/2 u a p 
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Therefore, we have in this subcase that 



(4.10) 1+ <<7sup 



\\d ai (aAc)\\ 



< c {e NM {t)} 1/2 v NAq (t). 



Secondly, we consider the region ©2 = {|£*| < |£| < !}• In this case, we can 
get that 



#122 <C / ^ 



ie - (,r ?0 («)M^fe)MT(e) ^xo^tf) 



Similarly as in f|4.9[) . 



ie - e*i 7 go(^)MT(^)M^(e) ^x^'KMO 



< c 



^X>2 



Hence, as in f|4. 10[) . the above two estimates imply 

#122 <C{£ NAq (t)} 1/2 V NJ , q {t). 
At last, we consider the region ©3 

= {M*\ < jlfUfl > !}• In this domain, > 

\\^\ and hence 

ie-e*r<cier<c(i + iei 2 + ie*i 2 )*. 

Thus, 



(4-11) #123 <C / dx |9 Q1 (a,6,c) 



^ + l£| 2 + l£*lT<7o(0) 



From Holder's inequality, 



< / (l + |£| 2 + |£*| 2 )^o(#)M^M^>^ M (£) a£XO 
(1 + |£| 2 + |^¥9o(0)M%)M%)wf /3M (O |^(0| 2 



< C 



l/1/2 W| ( 3|-£(0^ 2 2 'U2 2 • l/ 1/2 W| /3 |_f(C)wQ/3 
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where the computation similar to the estimate on the first factor in (|4.5[) by using 
|£| 2 + |£*| 2 = |£'| 2 + l£.| 2 instead of = has been performed. Therefore, 



I+<C I dx \d a ^a,b,c)\ 



<Csup ^ 2 w m _ e (Od^u 2 .\\d^(a,b,c)\\- v 1/2 w W -&)u aP 

< C{£ NAq (t)} 1/2 V NAq (t). 

Estimate on I 2 21- For the loss term, due to 



7 2 " 21 < J dx \d a *(a,b,c) 



we have by Holder's inequality that 

I 221 < C sup \d a2 (a, b,c 

X 

for \a 2 \ < N/2, and 

J 2,21 < C SU P 



v 1/2 u a p 



u 1 / 2 M q '/ 4 d^u 2 



\d a H^b,c)\\ 



v 1/2 u af3 



for |o;2 1 > N/2. Hence, in both cases, 

I 2ai <C{£ NAq (t)} 1/2 V NAq {t). 
For the gain term, by noticing that 

I+ 21 <C [dx \d a *{a,b,c)\ f |f-e.pffl)(0) 

■ M*(e*)M^(rHi^(e)|^> 2 (e:)^(o 

then, in the same way as for I 212 , when \a 2 \ < N/2, we have 
I+ 21 < C sup \d a *(a,b,c)\ • ^w^iOd^u, 

v 1/2 w\p\-e(t;)u a p < C {£N,i, q (t)} 1/2 V N . Lq {t). 



While when \a 2 \ > N/2, 



#91 < C ^ 



sup 
sup 



l|d Q2 (a,M 



v 1/2 u al 3 



x 11^(0,6,0)11 ■ ^ 2 w m _ e (0u a p 



over Di U 1 



over O3. 
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Therefore, for each c<2, 



lt,21<C{£N,L q {t)} 1/2 V NAq {t). 



Estimate on As for (|4.7|) . by [3T] Lemma 3, page 305], 

/2,2a <C{£ NAq (t)} 1/2 V NAq (t), 

and we also skip the proof of this term for brevity. 

Now, (|4.8j) follows by collecting all the above estimates. This completes the proof 
of Lemma 14.21 □ 



4.2. Estimate on £ • V x <pu. This subsection concerns the estimate on £ • V x (f>u. To 
this end, we first have the following result on the case of zeroth order derivative. 

Lemma 4.3. Assume that u and <p satisfy the second equation of (|2.3[) . It holds 
that 

(4.12) (h-V x cjm,u\ <^ t f R3 \ b \ 2 ( a + 2c ) dx 

+ C {\\(a,b,c)\\m + IIV^Hh! + HV^II • ||V X 6||} 
x \\\V x (a,b,c)\\ 2 + " 



x /2 { i_ p}m || 2 | 



c v 



(0 1/a {i-P}i 



Proof. Setting I3 = (^S, ■ V^u, u) and noticing u = Pu + {I — P}it = u\ 



u 2 , 



h = 



3 . m 



1 



1 



First, for /31, as in [12} 17], using the second equation of (|2.3|) to replace V x (f> yields 

h,i= [ V x (/>-b(a + 2c)dx = \— [ \b\ 2 (a + 2c)dx 
Jr3 2 at J M 3 

+ \b\ 2 1 1 V x • & + I • A(« 2 ) - 5 V K ■ & J 

+ / 6(0 + 2c) • {V. T (a + 2c) + V x • 0(u 2 ) - V^ajcfe, 

where 0(-) and A(-) are defined in (|2.2|) . Further by the Holder, Sobolcv and 
Cauchy-Schwarz inequalities, the above equation implies 



/3.1 < 



1 d 
2dt 



[ \b\ 2 (a + 2c)dx 
Jr 3 

+ C{\\(a, b, c)\\ H 2 + HV^II ■ ||V X 6||} ||V x (o, 6, c)|| 2 

+ C||(a,&, C )|| ff2 {||A( W2 )|| 2 + ||eM|| 2 } 



Noticing 



||A(n 2 )|| + ||eM||<C 
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then 13,1 is bounded by the right-hand term of (|4.12|) . Next, for 1^,1 and 13,3, one 
has 



ha < / |V X <; 



dx 



^- 1/2 wi 

<C [ |V»0| -\(a,b,c)\ 

<C||Vx^||hi f||V x (a,&,c)|| 2 



and 



/ 3 ,3 < C ||V x 0|| ioo (0 1/2 ^2 



(0 1/2 



"2 



Therefore, (|4.12j) follows from all the above estimates. This completes the proof of 
Lemma 14.31 □ 



For the case of higher order derivatives with respect to x variable, we have 
Lemma 4.4. Let N > 4, 1 < \a\ < N, and £>0. It holds that 



(4.13) 9" -e-V> 



<c||v^||^-i J2 {\\(0 1/2 ™-t, q (t,0d a {i-P}i 

1 <r l n I ^ at ^ 



l<|a|<JV 



\d a (a,b,c)\\' 



Proof. Let I4 be the left-hand term of (|4.13[) . Corresponding to u = Pu + {I — P}u, 

set 



(4.14) h=Y, 4m = ( d ° (U ■ V,<Mi) , wl^)d a u 



For i^i, one has 

4,1= E c£ v^- a v^i,^ M (o^ 

|qi|<|q| ^ ; 

< c E / |V.a Q - Q V|-|5 Ql («^,c)|-||9 Q u || i2 dx. 

Ql I < |0 

Here, when |cti | < N/2 and |ai| < |a|, the integral term in the summation above 
is bounded by 

snp\d^(a,b,c)\-\\W x d a - a ^\\-\\d a u\\, 

X 

and when \a\\ > N/2 or a\ — a. it is bounded by 

$vv\V x d a -^<t>\-\\d^{a,b,c)\\-\\d a u\\. 
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Thus, by the Sobolev inequality, £4,1 is bounded by the right-hand term of (|4.13JI . 
Similarly, we have for I^ 2 that 

h,2<C V f \V x d a - ai 4>\ ■ \\d ai u 2 \\ L , ■ \d a (a,b,c)\ dx 

\ ai \<\a\ jR3 & 

<c J2 su P ||a Qi u 2 || ir ||v E a Q - ai ^|| • \\d a ( ai b,c)\\ 

{\a 1 \<N/2}n{a 1 <a} X 

+ J2 8Up|Vxfl a - ai ^|-||fl ai «2|| •||fl a (o,6, C )||, 

{\u 1 \>N/2}U{a 1 =a} X 

which is also bounded by the right-hand term of (|4.13|) by Sobolcv's inequality. 
Finally, for ^4.3, since 



h.3<C V / \V x d a ^\- \^ 2 w. e (Od a -^u 2 



q 1 1 < a 



\e /2 ^i(Od a u 2 



dx, 



it can be further estimated as for I± t i and 14,2- Therefore, (|4.13[) follows from ()4.14j) 
by collecting all the above estimates. This completes the proof of Lemma T4. 41 □ 

For the case of higher order mixed derivatives with respect to both £ and x 
variables, we have 

Lemma 4.5. Let N > 4, 1 < \a\ + \j3\ < N with |/3| > I, and I > \/3\. It holds that 
1 



(4.15) ( d% \-i ■ - P}« , «ft Mi ,(f, - P}« 



\a\+\0\<N 



Proof. Due to 



where the summation is taken over a.\ + a 2 = a and /3\ + (3 2 = [3 with \[3i\ < 1, 
and for simplicity we denote each integration term in the summation as 1$. We 
now prove (|4.15[) by considering the following two cases: For the case |ai| < N/2, 
we have 

h<Csup\V x d a ^\- \t\ 1/2 w m _z{£)d%u 2 ■ \e /2 w\p\-i(Z)d%u 2 

x 

which is bounded by the right-hand term of (|4.15[) . For the case |ai| > N/2, it is 
easy to see that 



I5 < Csup 
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which is also bounded by the right-hand term of (|4.15[) . Therefore, (|4.5|) follows by 
combining both cases above. This completes the proof of Lemma l4~5l □ 

4.3. Estimate on V x (j> ■ V^u. Now we turn to estimate V x <f> • V^m. For the case 
with only derivatives with respect to x variable, we have 

Lemma 4.6. Assume -2 < 7 < 0. Let N > 4, 1 < \a\ < N, and £>0. It holds 
that 

(4.16) (d a (V^ • V*u) , w\ q (t, t)d a u) 

<C\\Vlct>\\ HN _A £ \\(0^\-^{I-P}u\\ 2 + 1^,(0,6,0)11^ 1. 

[|a|+|^|<iV,|/3|<l J 

Proof. Denote the left-hand side of (|4.16[) by Iq and write it as 

h = h,i + 2J C" ia2 {I 6t 2i{ai) + /6,22(ai)} 

ol\-\-ol-2— a, |ai I >1 

with 

/ 6 , 22 (ai) = (V x d ai cb- V^ a2u 2,wlM)d a u). 
We estimate term by term as follows. To estimate I§,\, notice 

+ (?) -27f e 2 9 (t)<C) 2 . 2g (t)V ? (C) 2 < City-^e 2 ^) 2 = C<f)«;2 / (f), 

where (£) > 1 and the fact that both g(i) = q + A/(l + t) and V^(£) are bounded 
by a constant independent of t and £ have been used. Then, from integration by 
part, 

<0 1/2 W£)d Q u || 2 , 

which is bounded by the right-hand side of (|4.16[) . For l6,2i((*i), it is straightfor- 
ward to estimate it by 

W«i) < C f \V x d ai <f>\ • |d Q2 (a,6,c)| • \\d a u\\ L2 dx 

< C||V^L w -i {||Vx(a,6,c)||^_ 1 + ||d Q u || 2 } . 
To estimate ^6, 22(^1), notice (£) 7+2 > 1 due to —2 < 7 < so that 



7 6>1 = (-gVx^- V C {^(C)},|5 Q M | 2 \ < Csup|V^|- 
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Thus, 

J 8 ,22(ai) < C [ \V x d ai <j>\ • ||(0^i-£(0V^ Q2 M2 || i2 • \\(0w-d0d a u\\ L2 dx 

<cpi<t>\\ HN _A J2 ||(o^i^i-^ 2 || 2 

[ |a| + |^|<JV,|/8|<l 

+ \\{0w-e(0d a u\\ 2 }, 

which is further bounded by the right-hand side of (|4.16l) . By collecting all the 
above estimates, it then completes the proof of Lemma T4. 61 □ 



For the case of higher order mixed derivatives with respect to the £ and x vari- 
ables, we have 

Lemma 4.7. Assume -2 < 7 < 0. Let I < \a\ + \(3\ < N with \/3\>l, and £ > \/3\. 
It holds that 

(4.17) (d$ {VA • V C {I - P}u) , tfp\- t , q {t, £) W - P}«) 

\a\ + \/3\<N 

Proof. Similar to the proof of Lemma 14.61 one can rewrite the left-hand side of 
KW, as 

(V^ • V^{I - P}u,wfa_ e (0d${I - P}u) 

+ E ^(v^fv^U-PK^i^i-pjtt). 

CK1+CK2— a,|cti|>l 

Then, similar to arguments used to deal with Iq i and ^6, 22(01) in Lemma 14.61 
(|4.17|) follows and the details are omitted for brevity. This completes the proof of 
Lemma 14771 □ 



5. A PRIORI ESTIMATES 

We are ready to deduce the uniform-in-timc a priori estimates on the solution 
to the VPB system. For this purpose, we define the time-weighted sup-energy 
X N ,i, q {t) by (dm]), and suppose that the Cauchy problem (fl75 |) -([T77 p of the VPB 
system admits a smooth solution u(t, x, £) over < t < T for < T < 00. We will 
deduce some energy type estimates on the basis of the following a priori assumption 

(5.1) sup X NAq (t) < 5, 

0<t<T 



where 5 > is a suitably chosen sufficiently small positive constant. 
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5.1. Dissipation of (a,b,c). Recall the fluid-type system (|2.3[) - (|2.4j) derived from 
the VPB system (|1.5|) - (|1.6j) and also (|2.2j) for the definition of the high-order mo- 
ment functions O(-) and A(-). In this subsection, we are concerned with the dissi- 
pation of the macroscopic component (a, b, c). 

Lemma 5.1. There is a temporal interactive functional Sffiit) such that 



(5.2) \£%\t)\ < Cl H 2 + Y, [\\d a {l-V}u{t)\ 

[ |a|<iV-l 

and 



\\d a W x (a,b,c)\[ 



(5.3) -£^{t) + Kl \\a\\ 2 + Y, \\d a V x {aAc)\? 

\a\<N-l 

<C Y p /2 <9 Q {I-P}t 

\a\<N 



+ C£ N j(t)l Y \\^ /2 d a {I-P}i 

[ \a\<N 

hold for any < t < T, where the simplified notion £N,i{t) = £N,e,q(t) with q = 
has been used. 



Y \\ d "( a ^' 

l<\a\<N 



Proof. Basing on the analysis of the macro fluid-type system (|2.3I) - (|2.4[) , the desired 
estimates follow by repeating the arguments used in the proof of [§J Theorem 5.2] 
for the hard-sphere case and hence details are omitted. Here, we only point out the 
representation of £ jy* (t) as 

Y I ^xd a c- k{d a {l-Y}u)dx 

\«\<N-l J ^ 

Y Y f ( d ^ ab i + d j9 a h - ?*«V„ • d a b) Qiji^il - P}u) dx 



\ a \< 

[ d a ad a \7 x -bdx 
\*\<n-i Jr3 

for a constant n > small enough. □ 

5.2. Construction of £N,t,q{t)- We are ready to prove the energy inequality (|1.19l) . 
In this subsection we consider the proof (|1.19j) in the following 

Lemma 5.2. Assume -2 < 7 < 0. Let N > 8, i > N and q(t) = q + -jj^jw > 
with 0<(}Cl,0<A<l and < $ < 1/4. Suppose that the a priori assumption 
(|5.ip holds for 5 > small enough. Then, there is £N,e,q(t) satisfying (|1.16|) such 
that 

d 



(5.4) 



dt 



£N,t,q{t) + KD Nt e tq (t) < 



holds for any < t < T, where T>N,i,q(t) is given by (|1.17j) . 
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Proof. It is divided by three steps as follows. Recall the VPB system (jl.5|l -( |1.6jl . 

Step 1. Energy estimates without any weight: First we consider the zero— th order 
energy type estimates. Multiplying (|1.5I) by u and integrating it over R 3 x R 3 gives 

\j t {\H 2 + llV^f} - (Lu,u) = (T(u,u),u) + V x cf>u,u\ . 

By applying (|1.14[) . (|4. 2[) and (|4. 1 2[) to estimate three inner product terms above 
and then using (|5.1|) . one has 



(5.5) ~ j|M| 2 + ||V^|| 2 - J \b\ 2 (a + 2c) rfxj + k \\^ 2 {I - P}, 

(0 1/2 {I-P}u 2 <C5V NAq (t). 



<C8\\y x (a,b,c)\\ z + 



CS 



Next we consider the energy type estimates containing only x derivatives. Applying 
d" with 1 < \a\ < TV to (jl.5p . multiplying it by d"u and then integrating it over 
R 3 x R 3 gives 

™ {\\d a u\\ 2 + l|a Q V^|| 2 } - (Ld a u,d a u) 



(d a T( Ul u), 8 a u) + ( d a [ -i ■ V x <jm ) , 8 a u ) + {-d a {V x ct> ■ V*u), d a u) 



After using (| 1 . 14[) . (|4.8p . (|4. 13[) . (|4. 16[) and then taking summation over 1 < \a\ < 
N , one has 



1 d 
Ydt 



Y {||a« u || 2 + ||9 a v^|| 2 } + K Y ||^ 1/2 a Q {i-Ph 



1<H<JV 1<|q|<AT 

< C\\(a,b,c)\\ HN ||V x (o,6 > c)||^_ 1 Y y /2 d a u\\+C{£* Aq (t)} 1/2 V NAq (t) 

l<\a\<N 



+ c ||v^|| HN _ 1 



\a\ + \P\<N,\P\<l 



+ \\V x (a,b,c)\\ H s-i >, 



which under the a priori assumption (|5.ip implies 
1 d 

l<|a|<JV " 1<|«|<JV 



Step 2. Energy estimates with the weight function w\p\—i iq (t, £): 
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Step 2.1. By applying {I — P} to the equation (jl.5|l . the time evolution of {I — P}u 
satisfies 

(5.7) d t {I - P}u + £ • V*{I - P}u + V x <t> ■ Vf{I - P}u + - P}u 

= K{1 - P}u + T(u, u) + ^- V x 0{I - P}u + [P, %\u, 

where [A £>] = AD — DA denotes the commutator of two operators A.D. and T,;, 
is given by 

T = i ■ V, + V x <j> • V e - ■ W x 4>. 

By multiplying (|5.7|) by w 2 _ l (t, £){I — P}u and taking integration over M 3 x M 3 , 
one has 

(5.8) ~||tu-/,,(t,0{i - PH 2 + <KO,^, g (*,OI{i - PW 2 ) 



1 d 2 



(t,0,|{i-P}t 



- p} u , wl Lq {t, e){i - p}«) + (r(«, «), ^(i, 0{i - P}«) 

+ (-V*0 ■ V e {I - P}u, ™ 2 f , 9 (t, fl{I - P}«) 



+ (-^-V^, W 2 (t,0|{I-Ph 



For the left-hand third term of ()5.8|) . notice 
1 d 2 - - A # 



(1 + t) 1 +* 



({P,T4u,w\ q (t,0{I-P}u). 



From (|2.ip . Holder's inequality in x integration and Cauchy-Schwarz's inequality, 
the right-hand first term of (|5.8[) is bounded by 



,1/2 



a,||xm<2c,(0" 7 '{i-p} 



z// 2 ™_ 



,(t,0{i-P}t 



< 2r/ 



^ 2 w^ q (t^){I~P}u 



C„ 



> x ' 2 {l-P}i 



for an arbitrary constant rj > 0. The right-hand second, third and fourth terms of 
(|5.8[) arc bounded by CST>N,e, q (t), where we have used (|4.3|) for the second term, 
velocity integration by part for the third term and also the a priori assumption 
(|5.1[) . For the right-hand last term of (|5.8|) . it is straightforward to estimate it by 



([p,r#,^(i,e){i-PH 



< rj 



- P}u|| 2 + C ri |||^ 1/2 V,{I - P}l 



\V x (a,b,c) 



cjv^i^J ^/ 2 {I-P} 



\V x (a,b,c) 
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for rj > 0. Plugging the above estimates into (|5.8[) and fixing a properly small 
constant 77 > yield 

2 



(5-9) ~\\w_ e ^0{l-P}u\\ 2 + K 



,1/2 



w_ A9 (t,0{i-P}« 



(1+*) 



c 



A/ 2 {I--p} u + v 1 / 2 V x {l-P}i 



\\V x (a,b,c) 



Step 2.2. For the weighted estimate on the terms containing only x derivatives, we 
directly use (| 1 .5[) . In fact, take 1 < \a\ < N, and by applying d" to (| 1 .5[) with 

Lit = L{I - P}w = -v{l - P}u + K{1 - P}u, 

multiplying it by w 2 _ e q (t,£,)d"u and integrating over M 3 x R 3 , one has 

2 



(5.10) l ~ \\v3- t ,&, t)d a 4 2 + (v(Od a {I - P}u, w\ q {t, Od a u) 

ld_ 

2 dt' 

= {Kd a {I - P}u, w\ q {t, £,)d a u) + (d a V x <}> ■ £,M 1/2 ,w 2 _ eq (t, £)d a i 

+ (d a r(u, u) - 9 Q (v^ • v e u) + d a • Vx^uJ , vLtJt, £)<9% 

For the left-hand terms of f|5 . 10[) . one has 

(v{0d a {\-P}u, w \ q {t,0d a u) 

v x > 2 W- t , q {t,C)d a {l -P}u 2 + (v(0d a {I - P}u,w\ q {t^)d a Pu) 



-w\ q {t,Z)Ad a u\ 2 



> 



and 



1 d 



2dt W -^^ d 



> 



a u\ 2 



Xi) 



lAPw-tafa Od a {l ~P}u -C \\d a (a, b, c)f 



'\(0w. e>q (t,0d a uf 



(1+t) 



l+i9 



2(1 



™ \\(t)w- t , q (t,t)d°{I P}u|r -C\\d a {aAc) 



For the right-hand first term of (|5.10p . one has 
(Kd a {I-P}u,w\ q (t,t)d a u) 



< 



(t, oo a {i - p}u\\ + c v \\x m <2c n (e) 7M) s a {i - P}«|| } 



v 1/2 w-e, q (t,0d a u 
V 1 / 2 w^ q (t,0d a {I-P}u 
+ C ri \jv 1 ' 2 d a {l-P} 



\d a (a,b,cW 
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where -q > 0. For the right-hand second term of ()5.10j) . we have that for 77 > 0, 



< Ci) 



v 1 / 2 d a u 
2 



c v ||a Q v a 



V*ff*{I - V}u\ + C v {\\a\\ 2 + ||V a (a,&,c)||^ N _i} 



For the right-hand third term of (|5.10[) , from Lemma 14.21 Lemma 14.41 and Lemma 
it is bounded by 



(5.11) CSV NAq (t) + 



cs 



(1+*) 



1+1? 



]T \\(Ow m _ itq (t,o{i-P}u\ 



*\ + \P\<N 



+ ||Vx(o,6,c)||^_ 1 } <C5T> N , e , q (t). 



Putting the above estimates into (|5.10|) . taking summation over 1 < \a\ < N and 
fixing a small constant 77 > give 



1 d 



l<\a\<N 



k Yl \\v 1/2 ™-i, q (t,0d a {I-Ph 

l<\a\<N 
K 



(1+*) 



1+1? 



\\{t) 1/2 w- t , q (t,t)d"t 



l<|a|<JV 



< C5V NJ . q (t) + C \ Yl \y /2 d a {l~P}u + \\V x {a,b,c)\\ 2 HN - 1 + \\a\f 

l<|a|<Ar 



Step 2.3. For the weighted estimate on the mixed x-£ derivatives, we use the 
equation (|5.7|) of {I — P}u. Let 1 < m < N. By applying dp with \/3\ — m and 
M + |/3| < N to (|5.7[) . multiplying it by g (t, ^)<9^{I — P}u and integrating 
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over R 3 xR 3 , one has 



1 d 



(5.13) l^-\\ Wm _e !q (t,Od%{I-F}u\ 



dp{ v {Od a {i - PM, wf 0l _ e Jt, omi - P} 



41™fW*>£)>IW- p }' 



dpKd<*{I - P}u, «;?.,_, ,(t, £)dpd a {I - P} 



r(«, u ) - v x <f> • v c {i - p}?i + ^ ■ v^{i - p}u 



+ ■ vj{i - p}«,^|_ li9 (u)9;{i - p>i 



For the left-hand second term of (|5.13j) . one has 
dp{v(Od a {I - P}u}, w 2 m _ e<q {t, Odpd a {I - P} 

\P1\<\P\ 



X\i\<2cMV m - l) d a {\-PY 



where rj > is arbitrary, and it is noticed that 



X|«|<2C,{C) 7(l ^- £) 9 a {I-P} W <C ^ 2 d a {I-P} 



For the right-hand first term of (|5.13j) . we have that for rj > 0, 



d Kd a {I - P}u, wf pl _ iA (t, £)d${I - P}u 



\Pi\<\P\ 



Crj 



Xm<2C n (0 im -' ) d a {l--P}r 



^ 2 w m ^ q (t,^{I-P}i 



<C V \\v 1/2 w\p\-t, g ( t >0dl{I-'Ph 

\h\<\P\ 



C v v l ' 2 d a {l-P}i 
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As in ([5. lip , from Lemma 221 Lemma T4. 51 and Lemma T4.7[ the right-hand second 
term of (|5.13|) is bounded by C6T>N,e,q(t). 
For the right-hand third term of (|5.13|) . 



^, e • vj{i - p} u , tf Phttq (t, - p}u 



E C Pl02 

<v 
+ c v 



-d^t ■ Vx%{I - P}u, w 2 m _ lA {t, (■)&${! - P}i 



E 

e| + |/3|<JV,|/3i| = |/3|-l 



v x/2 w lfi2l -t, g (t,Z)dl{I--p}u 



where rj > is arbitrary, and we have used 
w m-i, q (^0 = W\p\+i-e, q (t,0 w \f) 3 \+%-t,q( t ,Q 



The right-hand last term of (|5.13[) is estimated as 



-i,q\ 



(t,£)d£{I-P}n) < n v 1 ' 2 d a {l-P}i 



+ 0,(1 + 1^11^)^ E y /2 d a {l-P}u +||V a (a,6,c)||^-i 

for 77 > 0. Therefore, by plugging all the estimates above into (|5.13|) . taking sum- 
mation over {|/3| = m, \a\ + \/3\ < N} for each given 1 < m < N and then taking the 
proper linear combination of those N — 1 estimates with properly chosen constants 
C rn > (1 < m < N) and 77 > small enough, one has 



(5.14) 



1 d 

Ydt 



N 



J2C m J2 \\w m ^ q (t,0d^{l-P}u\ 



m=l \a\ + \P\<N ,\f3\. 



(l + t) 



E |^ 1/2 ^|_,, g (t,e)i^{i-p} 

\a\+\P\<N,\p\>l 

E ||«H^i-^(*»0^{i-p} 



1+1? 



\a\ + \P\<N,\P\>l 



<CSV NAq (t) + c\ Yl y /2 W- t , q (t,0d a {I-P}u +||V a (o,6 >C ) 

\a\<N 



\ H N-1 



Step 3. The proof of (|5.4[) : We take a proper linear combination of those estimates 
obtained in the previous two steps as follows. First of all, since 6 > is sufficiently 
small, the linear combination Mi x [ (|5.5[) + (|5.0[) ] + (|5.3p for Mi > large enough 
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gives 



(5.15) 



d I M x 

dt 1 T~ 



V (\\d a u\\ 2 + \\d a V x <f>f)- [ \b\ 2 (a + 2c)dx 



\a\<N 



; Y /2 d a {l-V}u +\\a\\ 2 +\\V x (a,b,c)\\ 2 HN -A <C5V NAq (t), 

\a\<N J 

where the corresponding energy functional is equivalent to IMI^^n) + || V^H^n . 
Now, the further linear combination M 3 x [M 2 x (|STT5| + ([S3)) + ([532]) ] + (|5~TI)) for 
M2 > and M3 > large enough in turn gives 

^NJ,q(t) + KV Nt £ iq (t) < C5V NAq (t), 

where £jy,^ g (t) is given by 
(5.16) £ N ,i, q (t) 

= M 3 





Mi 


M 2 < 






{ 2 



E 

\a\<N 



u|| 2 + II^V^^H 2 ) - / \b\ 2 (a + 2c)dx 



l<|a|<JV 
TV 

+E c ™ e ihi/ji-/,,(*,OT{i-p}«ir- 

"i=l |a| + |/9|<JV,|/8|=m 

Notice (|5.2I) . It is easy to see that 

ZNAlit) ~ H/V,£,g(*)- 

Therefore, f)5 -4[) follows from the above inequality since J > is small enough. This 
completes the proof of Lemma 15.21 □ 



6. Global existence 



Recall (|1.18[) for XN,e, q (t). To close the energy estimates under the a priori 
assumption (|5.1[) , one has to obtain the time-decay of £N,e-i,q(t) and || V^H^-i ■ 
The key is to prove the following 

Lemma 6.1. Assume —2 < 7 < and J R3 clq(x) dx = 0. Fix parameters N , £q, t 
and q(t) = q+ (i+ t \4 as stated in Theorem ll.il Suppose that the a priori assumption 
()5.1|) holds true for 5 > small enough. Then, one has 

XN,e,q(t) < C {^N.e.q + ^N,i,q(t) } , 



(6.1) 



for any t > 0, where eN.e.q is defined by 



(6.2) e NAq = J2 \\ w W-ui ^)9pU \\ + (1 + N + {0~^>o 

\a\ + \P\<N 



VLASOV-POISSON-BOLTZMANN SYSTEM 



As preparation, we need an additional lemma concerning the time-decay esti- 
mates on the macroscopic quantity (a, b, c) in terms of initial data and XN,t(t). 
Here and afterwards, for notational simplicity, we write X^^(t) = X^^^ q {t) when 
<7 = 0, and in a similar way, cn,i = £N,e,o and £N,i{t) = £N,£,o(t) will be used. 

Lemma 6.2. Under the assumptions of Lemma \6.1[ one has 

(6.3) ||u(t)|| + ||V X 0(*)|| < C(l + t)-i {e N/ + X N/ (t)} , 
and 

(6.4) || V x w(t) + ||V^(i)|| <C(l + t)-i{e N ,i + X Nti (t)}, 
for any < t < T . 

Proof. By Duhamel's principle, the solution u to the Cauchy problem (|1.5p - (|1.7p of 
the nonlinear VPB system can be written as the mild form 



u(t) = e tB ito 



,(t-s)B 



G{s) ds 



with G = §£ ■ V x <p u — V x • V^u + T(u, u). Notice that PoG(i) = for all t > 
and the condition (|3.3p holds. Applying Theorem 13. 11 one has 



(6.5) Kt)|| + ||Vx0(*)|| 

<C(l + t)-i( (0~ 



and 



>U 



+ C + 



■||(l + N)ao|Ui 



ds 



•6) ||V a u(t)|| i2firi 



|V^(*)|| 



<C(l + i)-4 <0~*'°«o 



(0-^°G( S ) 



if(Hi'- 2 ) 



+ ||(l + |x|)o ||j 



(0^ 4 V X G( S ) 



ds, 



where we recall that £q > 2<j\ = 2(| + |) = 5/2 is a constant. In what follows, we 
shall estimate the right-hand time integrals in the above two inequalities. 
To do that, we claim that 



(6.7) 



£ \kr iea d a G(t) 



< C£ 



(*), 



|a|<JV-l 



for any < t < T. Here, recall £ > 1 + max|vV, 4^ — i|. In fact, one has the 
following estimates. 
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Estimate on terms containing 4>: Direct calculations yield 



(6.8) 



(0- 



-V x (f> ■ Vju + ^£ ■ V x (f>u 



< c 



(0- 



|v x ^||ij{||v € «|Uj + <oil«IUs} 



< 



G]|Vx^||{|^}-^°v^ 



<0" 



} 



< c 



c||v x 0|| {|w-|A,(t,0Vc«| + ||t«_^ + i(t,0«||} 



Here, by recalling the assumption I — 1 > % — A, (|6.8|l is further bounded by 
C£N,i-\\t)- 

Similarly, by applying L^-norm to the low-order derivative terms and using 
Sobolcv's inequality, one also has 



j2 (^-i^j-V^-Vcu+^V^u] 

\<N-1 



<C\\W x <f>\\ H N hm-v-vMVpA <CE N ,e-i{t)- 

\a\ + \0\<N 

Estimate on terms containing T(u, u): First consider 
(6.9) (£>-i*T(«,tO v ^C^hr^ ^ (\\u\\ Ll ,\\u\\ Ll ) 

where Holder's inequality has been used. For the loss term, by noticing 

T- (\\u\\ Ll ,\\u\\ L 2) = \\u(0\\l* jj |£-£*P?o(0)M%)|K£OIU^*^ 

<C|K0IU r (C) 7 sup|| U (e*)|U Sj 



\u\\ Ll ,\\u\\ Ll ) <Csup|Ke*)IUr (0 7( ^ +1) ? 



L- 



one has 
(6.10) 

<C(||V f «|| + ||V|«||) w_^ +1 (t,0u <C8 NA ^{t). 

For the gain term, 

JiM ■■= (0-^r+(||«|| £ ;,||u|| L j) 

= @~ ito JJ ie-e*i 7 9oWM^(^)ime:)ii i2 j| W (e')iiL^e*^. 

We then consider the following three cases: 
In Pi={|g.|>j|g|} : 

Mt,o < [ ie-€.r<ft,(fl)M^«.)M^(on«(e:)iu»n««')ik- 
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From Holder's inequality, as before, Ji(i, £) is bounded by 

c^-^+hU ie-^rffl,((?)M^(^)M^(oiiu(e:)iii;iiu(oiii; 

so that 

l|Ji(t)x Dl ||| ? < c / ie-^l 7 goWM^(^)M^(oll^(e)HMk(Oll!= 

<C [ \^-^PM^(^)M^mu(^)\\hMO\\h 

< Csu P |K£*)ll 2 L2 • / \u(x,o\ 2 < ce N ^-i{tf 

£. x 

In B 2 = < |£| < 1}: From Holder's inequality, one has 

Ji(U)xo 2 < { / x^-^qo{e)M^{^)\\u{0\\l l \Hmii 

\ J £,UJ 

In ©2, |£ — > an( i moreover, 

i€ , i<2i^i + i^i<|iei ) iei<2i^i + iei<2i$i. 

These imply that in ©2, 

i£-&r <cmin{ie'r,iei 7 }xk'i<3,iei< 3 - 

Then, using the above estimate and then taking change of variable (£, £*) — > (£', 
one has 

\\Ji(t) X ®J%<C [ X a ^-t*rqo(0)MH^)\\u(Q\\l4u(e)\\h 

<c / xici<3,ic.i<3nun{|cr,i^r}ii««.)iii»ii«(oiii ; 

In D 3 = < |£| > 1}: In this case, as in the proof (|4.11|) for 7^123' f rom 



Holder's inequality, we can deduce that J\{t, £) is bounded by 

c<£}-^° ( / (i + iei 2 + ikiTsoWM^iKOiiLiwoii? 



which further implies that ||Ji(t)XD 3 \\ 2 L 2 is bounded by 

c f (0- 7i ° (i + l£l 2 + l£*l 2 ) 7 gb(0)M*(&)NOIl£ a IKOI& 

•/£,£,,w 

< c / {(e'}-^ + <£)- 7 '°} (i + l£l 2 + l£*l 2 ) 7 9 oWm%)hoh| 2 KOliij 

< c / + <e*r 7 '°} (i + \e + ie*i 2 ) 7 n«(6)iiij ii«(onis- 

J£,£. 
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Thus, it follows that 

\\Mt)x 03 \\l t <cj (o-^inmh/t- J iiu&)iiij«. 

+c|«.>-' r * ,+7 Kwiii s de.- / hoii!^ 

<C w_ia +i (t,Ou 2 IKII 2 < C£ N ,t-x(ty 

2^2 

Collecting the above estimates, one has the estimate on the gain term 

Combining the above with (|6. 10[) and recalling (|6.9[) . we have 
(6.11) 

Now, we consider the L -norm of {i)-i la d a T(u,u) with \a\ < N — 1. Notice 



(0-^°r( u , w ) <C£ Nlt (t). 



hr ¥o d a T(u,u)\\ <c J2 E E||^)"^ or± ( 9aiu ' 5a2M ) 



|«|<JV-1 



\a\<N-l Qi+Q2=a ± 



|q|<W-1 ai+a 2 =a,|ai|<A'/2 ± 

+ c E E El^-^ridi^uiuj.iiv^uiiHi; 

|q|<JV — 1 Qi+ct 2 =a,|ai |>7V/2 ± 

where Sobolev's inequality < C|| V^/H/fi has been used. Then, as for (|6.9p . 

one has 

E ||(0"^°9 a r( u , u )| < c?£^-i(t). 

|q|<A-1 

and (jlTTTj) give flB7T|). □ 



This together with 

Hence, (|6 . Tf) together with (|1.18[) imply 



TV, £-1 



\a\<N-l 



< C(l + S )~ sup (1 + fl)»£j Vl /_l(«) < C(l + fl)-5XjV,i(t), 
0<s<t 

for any < s < t < T. Further plugging the above estimate into (|6.5[) and (16. 6[) 
yields (|6.3[) and (|6.4p . respectively, where the following inequalities 

ft 

(l+i-s) _ z(l + s)-3ds < C(l + t)~*, 



(i + t-sy*(i + s)-5ds < c(i + ty*. 



together with (|1.16l) and (|5 . for the definition of £-N,l,q have been used. This 
completes the proof of Lemma 16.21 □ 
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Proof of Lemma l6.lt We divide it by three steps. Recall (|1.18[) for XN,e,q(t). 
Step 1. From Lemma l5~2l the time integration of (|5.4|) directly gives 



3.12) 



sup £ N ,e, q (s) < £ N ,i, q (0) < Ce 2 N t q 



Q<s<t 



Here, noticing A x (j)(0,x) = ao(x), we have removed || V x (j>(0, x)\\^ N from £/v,^ 9 (0) 
by the definition of (jl.lpp or (|1.16|l in terms of the inequalities 



and 



|V^o|| <C\\a \\% 3 \\a \\% 3 



^o|| ffN -i = || 1 a || JfJV _ 1 < C||ao||ifj> 



Therefore, £N,e, q (0) < Ce 2 N e holds true with CN,e,q given by (|6.2j) . 

Step 2. Take < e < 1/2 small enough. Recall (fOj) . Notice that (f5~4"l) also holds 
true when I is replaced by I — 1 since all the conditions of Lemma 15.21 are still 
satisfied under the assumption that I > 1 + N and sup 0<s<T X/v^ ig (s) < 5 with 
5 > small enough. Thus, it holds that 

-jj£N,l-l,q{t) + KV N ^-l,q(t) < 0. 

Multiplying the above inequality by (1 + i) 3 / 2+e gives 



(6.13) 



dt 



(1 + t)* +e £ N: i-i, q {t) + k(1 + t)i +t V Ni e- h q{t) 



< ( 2+ e )( 1 + f ) 



'N,t-l,q 



,*(*)• 



Similarly, from (|5.4[) with £ replaced by I — 1/2 and further multiplying it by 

(1 + i) 1 / 2+e , one has 



3.14) 



d r 

dt 



(1 + t)~^£ N/ _ hq (t)\ + k(1 + t)-^V N4 _ hq (t) 



Observe from (fTTU|) . ([Llo) and (fTTT7|) that 



W*) + ii(M(*)ii + nv^wir > ^_i, 9 (t) 



holds for any given ^. Then, from taking the time integration over [0,t] of (|6 . 1 3[) . 
(I6.14p and (|5.4[) and further taking the appropriate linear combination, one has 

(1 + t)i +t £ N ^ q {t) < C£ N ,i. q {0) 

+ C /"*(! + =s)i+ e {||(fo, c)( S )|| 2 + HV^^Cs)!! 2 } ds. 



Here, applying the first estimate (|6.3p in Lemma 16721 to the right-hand time integral 
term of the above inequality and noticing 

(i + s )i +e (i + ds < c(i + ty, 
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it follows that 

(1 + t)-2+'£ N j^ q (t) < CS NAq (0) + (7(1 + ty {e 2 N j + X N , e (t) 2 } , 
which implies 

(6.15) sup {l + s)iS N ,e-i, q (s) <C{e 2 Ne „ + X N . e (t) 2 } 

0<s<t ' ' 

holds for any < t < T. Here, £N,e,q(0) < Ce 2 N l q has been used. 
Step 3. Notice 

llV^H^-x < (7{]|V^|| 2 + ||V,a||^_ 2 } < cjllV^f + ||V^||^ (ff „_ 2) | . 
Then, it is direct to deduce from (|6.4j) 

(6.16) sup {l + s)i\\Vl<l>{s)f HN -, < C{e 2 Ne + X NJ (t) 2 } . 

Q<s<t 

Hence, the desired estimate (|6.ip follows by summing (|6.12p . f|6 . 1 5[) and (|6 . . 
This completes the proof of Lemma 16.11 □ 

Now, we are in a position to give the 

Proof of Theorem ll.lt Given -2 < 7 < 0, we fix N, £ and q(t) = q + jj^w as 
stated in Theorem ll.il The local existence and uniqueness of the solution u(t, x, £) 
to the Cauchy problem (|1.5|l - (|1.7|l can be proved in terms of the energy functional 
£N,e,q(t) given by (|1.16[) . and the details are omitted for simplicity; see [35]. Then, 
one only has to obtain the uniform-in-time estimate over < t < T with < T < 
00. In fact, by the continuity argument, Lemma |6 . 1 1 implies that under the a priori 
assumption (|5.ip for 6 > small enough, one has 

(6.17) X NAq (t) < Ce% Aq , < t < T, 

provided that £N,l,q defined by (|6.2|) is sufficiently small. Recalling the condition 
(|1.12p for initial data Ug which coincides with (|6.2p . the a priori assumption (|5.ip 
can be closed. Then, the global existence follows, and (| 1 . 1 3|) holds true from (|6.17D . 
([TT51) . HHg) and (fTTTDj) . This completes the proof of Theorem O □ 
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